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Abstract 


We apply the canonical perturbation theory to the semi-quantal hamil- 
tonian of the SU(3) shell model. Then, we use the Einstein-Brillowin-Keller 
quantization rule to obtain an analytical semi-quantal formula for the en¬ 
ergy levels, which is the usual semi-classical one plus quantum corrections. 
Finally, a test on the numerical accuracy of the semiclassical approximation 
and of its quantum corrections is performed. 


PACS: 03.65.Sq; 05.45.+b 
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In the last few years, there has been considerable renewed interest in the 
semi-classical approximation, due to the close connection to the problem of 
the so-called quantum chaos [1,2]. One important aspect is the semi-classical 
quantization formula of the energy levels for quasi-integrable systems [3,4], 

It has recently been shown [5,6] that, for perturbed non-resonant har¬ 
monic oscillators, the algorithm of classical perturbation theory may also be 
used in the quantum-mechanical perturbation theory, with quantum correc¬ 
tions in powers of h. 

In this paper, on the contrary, we calculate the quantum corrections to 
the semi-classical quantization [3,4] of a many-body model related to nuclear 
physics. Its classical counterpart, obtained in the limit of the number of 
particles that goes to inhnity, is represented by a non-integrable hamiltonian 
with two degrees of freedom [7,8,9,10]. The semi-classical quantization of this 
model has been studied in [7] and here we calculate the quantum corrections 
and then analyze their numerical accuracy. 

The model is a three-level schematic nuclear shell model, whose hamil¬ 
tonian is: 

H = 't + Gl„ ( 1 ) 

k=0 ^ fc^z =0 


where 


M 

Gkl = ^ 

m=l 


( 2 ) 


are the generators of the SU(3) group. This model describes M identical 
particles in three, M-fold degenerate, single particle levels e*. There is a 
vanishing interaction for particles in the same level and an equal interaction 
V for particles in different levels. We assume 62 = —eo = e = 1, ei = 0. 
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For the SU(3) model the semi-quantal hamiltonian [11] is dehned as [8]: 

H{pi,P 2 ,qi,q 2 ]M) =< qiPi, q 2 P 2 ] M\ — \qipi, q 2 P 2 , M >, (3) 

where \qiPi, q 2 P 2 ] M > is the coherent state, given by: 

\qiPi, q2P2] M >= exp [ziGm + 2 : 2 ^ 02 ]|00 >, (4) 


with: 


Mk + iPk) = 


^l + zlz, + zW ^ 

and |00 >= n^^aQ^|0 > is the ground state. Here 1/M plays the role of the 
Planck constant h [10]. 

As discussed in great detail in [10], the semi-quantal hamiltonian is: 
H{pi,P 2 ,qi,q 2 -,M) = -1 + ]^{p\ + ql) + (^2 + ^ 2 ) + ^x[l “ 


X[(91+92) -{P1+P2) -( 9 i-Pi)( 92 -P 2 )- 49 i 92 PiP 2 - 2 (gi+g 2 -Pi-P 2 )], ( 6 ) 


with X = MV /e. The phase space has been scaled to give {qi + q^+Pi+pi) + 
2. The classical hamiltonian can be obtained in the ’’thermodynamical” limit 
[ 10 , 12 ]: 

^d(Pi,P2,9i,92) = lim M(pi,P2,9i,92;M), (7) 

M—^OO 

and the semi-quantal hamiltonian is given by: 


H{pi,P2,qi,q2]M) = Mci(pi,p2, 9i, 92) + ^fgc(Pi,P 2 , 9i, 92; M), (8) 

where Hqc is the hamiltonian of quantum corrections. 

Through the canonical transformation in action-angle variables [11]: 

qk = \/^k cos {6k), pk = sin (Ok), k = 1,2 (9) 
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the semi-quantal hamiltonian can be written: 


h, 01,02; M) = Hoih, h) + xV{h, h, 01, 02; M), (10) 

where: 

Ho{h,l 2 ) = -l+Ii + 2l2, (11) 

V{h, h, 01, 02; M) = cos (20 i)+/2 cos (202)]+/i/2 cos (202 - 20i). 

( 12 ) 

We applied a canonical transformation {Ii, 12 , 01 , 62 ) {Ii, 12 , 61 , 92 ) in 
order to obtain a new hamiltonian that depends only on the new action 
variables up to the second order in a power series of y: 

H{h, h; M) = Hoih, h) + xHiih, h M) + x^H 2 {h, h M). (13) 

It is well known that the canonical perturbation theory presents many dif- 
hculties which are essentially related to the so-called small denominators. 

The resonance of the unperturbed frequencies uji = = l,a ;2 = = 2: 

mooi + nuj 2 = 0, (14) 

can lead to divergent expressions in the perturbative solution to the prob¬ 
lem. This drawback occurs only if the integer numbers m and n are present 
as Fourier harmonics in the perturbation theory. We will show that the 
resonance condition (14) is not satished up to the second order in y. 

We assume that the generator S of the canonical transformation may be 
expanded as a power series in y: 

S{ii, I 2 , 01, 02; M) = /i0i-|-/202-I-x5'i(/i, I 2 , 01, 62 ] M)+x‘^^ 2 ( 11 , 12 , 0i, 62 ] M). 

(15) 
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2882 

(16) 

06k 

^k + X 

oik 

2^82 

(17) 

+ ’ 
oik 


The generator S satisfies the equations: 

dS 

~ _ dS 

Uh — —^ 

dh 

with /c = 1, 2. From the Hamilton-Jacobi equation: 


dS dS dS dS ~~~ - 

^ Hq{Ii, l2)+Hi{Ii, /2; M)+Fr2(/i, h', M), 

(18) 

we have a number of differential equations obtained by equating the coeffi¬ 
cients of the powers of y: 


dSi dSi' 


dOi 


dOo 


1 + A + 2 / 2 , 

(19) 

+ V{ii,i2,ei,e2-,M), 

(20) 

fdVdSi dVdSA 
U/i dOi + dh 002 I 

(21) 


The unknown functions Hi, Si, H 2 and S 2 may be determined by averaging 
the time variation of the unperturbed motion. At the first order in y we 
obtain: 


~ ~ ~ 1 r 27 r r27r 

Hi{ii,i2-,M) = — / deide2V{ii,i2,ei,e2;M) = o, (22) 

4:71^ Jo Jo 


and 


Si{ii,i2,ei,e2-,M) = - ( 23 ) 

(Kn)} imu;i + nu;2) 
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where {(m, n)} = {(2, 0), (0, 2), (—2, 2)} are the Fourier harmonics of the 
perturbation potential V. The resonance condition is not satished, and we 
have: 


= -^[1 - - / 2 )(/isin( 20 i) + 

2 ^^ ~ ]^]-^ 2 sin ( 26 ' 2 ))] - -[1 - —]/i /2 sin ( 26^2 -26'i). 


At the second order in 


^2(/i,/2;M) = 


/*27r /*27r 


dTT^ 


'0 JO 


d9id92 


'dVdSi dVdS: 


+ 


dll ( 96*1 < 9/1 ( 96*2 


— ^[1 ~ + -^1 + 2/2)(2/1 — All + -^2 — iih — li , 


and: 


(24) 


(25) 


^2(/i,/2,^1,^2;M) = - ^ ^ ^ (26) 

{im,n)} irnu;i + nuj2) 

where {{m, n)} = {(2, 0), (4, 0), (2, -4), (4, -4), (2, -2), (0,4), (2, 2)} are the 
Fourier harmonics of the function IF, given by: 


fF(Ji, J 2 , 9i, 92; M) = H2{h, h; M) - 


'dVdSi dV dSi 

+ 


dll ( 96*1 ( 9/2 ( 96*2 / 

In this case too, the resonance condition is not satisfied and we have: 


(27) 


*S' 2 (/i, - 62 , ^ 1 , ^ 2 ; Af) — -[1 —[3/i/ 2(1 — /i — I 2 ) sin (26*i) + 

o IVl 

+Ji(l — 3/1 + 2/^ — 2/2 + 3 / 1/2 + I 2 ) ( 4 ^ 1 ) + 

+3/i/2(/i + li - 1) sin {29i - 4 ^ 2 ) + 

+/i/ 2 (/ 2 -/i)sin (401-402) + 
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+3/i/ 2(1 — /i — I 2 ) sin {26i — 2 ^ 2 ) + 

+—/ 2(1 — 2Ji -\- Ii — 3/2 + 3 / 1/2 + 2 / 2 ) sin ( 4 ^ 2 ) + 

+ 3 / 1 / 2 (/i 4“ -^2 — 1) sin ( 2^1 + 2 ^ 2 )] (28) 

In conclusion: 

H(h,h\M) = -lVh+2h+^[l-^\(-l+h+2U)(2h-iil+'h-hh-il)^ 

(29) 

This approximate semi-quantal hamiltonian depends only on the actions. 
Thus, a semi-quantal quantization formula may be obtained by applying the 
Einstein-Brillowin-Keller rule [2,3]: 

/fc = + 2 ) /c=l,2 (30) 

where 1/M plays the role of the Planck constant h. In this way we have: 


= EZJM) + K.nSM) 


(31) 


where: 

K\n,{M) = -l+(ni + ^)^+2(n2+^)^+^[-l+(ni + ^)^+2(n2+^; 


1. 1 


132 1 


1. 1 


1. 1 




M2^’ 

(32) 


is the semi-classical quantization formula, and 


Kln.iM) - -^[-1 + (ni + -)— + 2(n2 + -)—] 


1 . 1 


1 , 1 


X 


4M^ ' 2'M ' ^ ^ ' 2'M- 

xl2(n.+i)i-4(n,+i)=jL+(„,+ l)^_(,„.+ l)(„,+ l)^_(„,+ l): 


M2J’ 

(33) 



are the quantum corrections. 

In order to test the accuracy of the semiclassical approximation and its 
quantum corrections, the eigenvalues of the hamiltonian (1) must be calcu¬ 
lated. A natural basis can be written: \bc >, meaning b particles in the 
second level, c in the third and, of course, M — b — c in the hrst level. In this 
way 100 > is the ground state with all the particles in the lowest level [7,8]. 
We can write the general basis state: 

\l>c>= Ai\00>. (34) 

where is the normalizing constant. 

We can calculate the expectation values of ^ and, therefore, the eigen¬ 
values and eigenstates of In this way, the energy spectrum range is 
independent of the number of the particles: 

< b'c\^\bc>= ^(-M + 6 + 2c)(5f,y(5,y (35) 

where: 

Qyc',,c = i/Kb-l){M-b-c+l)iM-b-c + 2) 

+ \/(b + l){b + 2){M - b - c){M - b - c - l)5b+2,6"^cc' 

+ \J c(c — 1)(M — b — c + 1)(M — b — c + 2)(5ft 
+\/ic+l)ic + 2){M -b-c){M-b-c- l)5fe,y5c+2,c' 

+ ^J{b -h 1)(6 -h 2)c(c - l)5fe+2,fe' bc- 2 ,c 
+ \/b{b - l)(c -h l)(c -h 2)(5^_2^y(5c+2,c' (36) 
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and X = MV/e. The expectation value < > is real and symmetric. 

For a given number of particles M, we can set up the complete basis state, 
write down the matrix elements of < ^ > and then diagonalize < ^ > 
to hnd its eigenvalues. < ^ > connects only states with A6 = —2,0,2 
and Ac = —2, 0, 2 which makes the problem easier. We group states with 
b,c even; b,c odd; b even and c odd; b odd and c even. This means that 
< ^ > becomes block diagonal containing 4 blocks which can be diagonalized 
separately. These matrices are referred to as ee, oo, oe and eo (for further 
details see also [17]). 

Then we compare these ’’exact” levels to those obtained by the semi- 
quantal perturbation theory. A very good agreement is displayed (see Fig. 
!)• 

In Table 1, we show the difference between the ’’exact” levels and those 
obtained by the semi-classical and semi-quantal perturbation theory. We 
observe that the algorithm provided by the semi-quantal perturbation theory 
gives better results than that of the ordinary semi-classical perturbation 
theory. 

Obviously if 1/M, no matter how small, is kept hxed, this semi-quantal 
approximation on the individual levels has the meaning of a perturbation 
theory in 1/M [5,6,13]. Therefore, the accuracy of the approximation de¬ 
creases for higher levels [14]. To obtain a better agreement it is necessary, as 
is well known, to implement the classical limit 1/M—>0, s^cx) and, at 

the same time, to keep the action = {uk + l/2)/M constant [15,16]. 

Finally, we stress that, for systems with a hnite number of Fourier har- 
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monies, like the SU(3) model, rational frequencies do not give rise to the 
problem of small denominators up to a certain order of the canonical pertur¬ 
bation theory. 
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Table Captions 


Table 1: The differences for the first 10 levels, with x : 
for the eo class, are the ’’exact” levels, are the 
and are the semi-quantal levels. 


= 0.75 and M = 100, 
semi-classical levels. 
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Figure Captions 


Figure 1: Comparison between ’’exact” levels (left) and those obtained by 
the semi-quantal perturbation theory (right); with y = 0.75 and M = 100 
for the eo class. 
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j^ex _ j^sc\ \j^ex _ j^sq 


1.3086796-10-3 

1.6146898-10-3 

2.7745962-10-^ 

1.4380813-10-3 

1.5463829-10-3 

1.0503531-10-3 

1.9035935-10-3 

4.4906139-10-^ 

6.0987473-10-^ 

1.8021464-10-3 


1.1860132-10-3 

1.3964772-10-3 

2.2178888-10-^ 

1.1500716-10-3 

1.3815761-10-3 

7.1579218-10-^ 

1.6558170-10-3 

3.5130978-10-^ 

2.4980307-10-^ 

1.4950633-10-3 


Table 1 


17 




